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ABSTRACT 

Despite  gasdyneunic  non-idealities  in  the  flow  pro- 
duced in  a shock  tube,  pressure  measurements  at  three 
different  locations  in  the  driver  section  of  the  shock 
tube  revealed  that  the  expansion  wave  generated  in  rel- 
atively weedc  expansions  could  be  viewed  effectively  as 
a simple  centered  expeuision  fan  after  am  empirical  shift 
of  the  actual  origin  of  the  expansion  wave  to  a "virtual" 
origin.  The  resulting  centered  expansion  fan  was  used 
to  study  at  two  locations  the  condensation  of  H2O  and 
D2O  vapors  in  an  excess  of  the  carrier  gas  argon,  with 
simultaneous  pressure  amd  light  scattering  measurements. 
The  isentropic  flow  within  the  centered  expansion  fan 
was  found  to  be  preserved  up  to  the  point  of  the  detect- 
able onset  of  condensation  by  tailoring  the  onset  con- 
ditions to  occur  at  the  tail  of  the  expansion  fan,  thus 
rendering  a simple  analysis  of  the  experiments  possible. 
The  onset  conditions  for  H2O  vapor  were  found  to  be  in 
agreement  with  previous  findings  in  supersonic  nozzles 
and  shock  tubes,  amd  they  were  well  predicted  by  the  so- 
called  classical  theory  of  homogeneous  nucleation.  The 
condensation  of  D2O  vapor  was  found  to  exhibit  similar 
trends  as  those  of  H2O  vapor  condensation  despite  the 
slight  differences  in  physical  properties  between  them 
due  to  isotopy. 
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ct  classical  theory  of  homogeneous  nucleation 
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INTRODUCTION 


In  addition  to  the  well-known  cloud  cheunber  and 
supersonic  nozzle  methods  of  condensation  research,  the 
shock  tube  offers  a technique  that  has  many  desirable 
features.  As  a batch  experiment  like  the  cloud  chamber 
it  has  the  advantage  of  accurately  controlled  experi- 
mental conditions  amd  easy  exchange  of  test  substances. 
Only  small  amounts  of  condensing  substances  are  required 
which  can  be  prepared  with  great  care.  Thus,  toxic  as 
well  as  costly  materials  can  be  studied.  The  wholly 
confined  system  can  be  thoroughly  cleaned  by  pumping  it 
to  high  vacuum  prior  to  an  experiment.  This  greatly  re- 
duces the  problems  of  contamination  and  seeding  by 
foreign  substances  which  may  lead  to  undesirable  compli- 
cations such  as  binary  nucleation  or  heterogeneous  nu- 
cleation,  particularly  at  low  cooling  rates.  With  appro- 
priate cooling  and  heating  facilities  a wide  range  of 
substances  may  be  investigated  over  a wide  range  of 
experimental  conditions. 

The  shock  tube  was  first  applied  to  condensation 
studies  by  Wegener  and  Lundquist* ; and  the  first  streaJc 
photographs  showing  condensation  zones  in  the  expansion 
fan  in  the  driver  section  were  provided  by  Glass  and 
Patterson* . Since  then  there  have  been  many  different 
applications  of  different  parts  of  the  shock  tube  flow  to 
condensation  studies.  For  exeunple,  Homer*  »**  utilized  the 
constant  thermodynamic  condition  behind  the  incident  shock 
wave  to  study  the  nucleation  and  growth  of  lead  particles. 
Kung  and  Bauer*  investigated  the  nucleation  of  iron  vapor 
with  the  so-called  tailored  interface  technique  which 
exploited  the  expansion  fan  produced  as  a result  of  the 
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tailored  interaction  between  the  reflected  shock  and  the 
contact  surface.  The  method  described  in  this  work  makes 
use  of  the  cooling  capabilities  of  the  unsteady  isentropic 
expansion  wave  in  the  driver  section  to  induce  condensa- 
tion. The  onset  of  condensation  as  indicated  by  the  first 
appearance  of  condensation  aerosols  can  be  detected  by 
Rayleigh  scattering  of  laser  light  whose  wavelength 
(X  = 4416  % for  a Helium- Cadmium  laser)  is  much  larger 
than  the  size  of  the  particles  generated  under  these  experi- 
mental conditions  (normally  a few  hundred  Angstroms);  while 
the  sudden  release  of  the  latent  heat  of  condensation  gives 
rise  to  a pressure  "bump".  Such  a method  has  been  used  by 
Barschdorff*  to  study  the  carrier  gas  effects  on  the  homo- 
geneous nucleation  of  water  vapor.  Courtney’  has  also 
used  the  expansion  wave  to  study  the  condensation  of  water 
vapor,  methanol  and  carbon  tetrachloride.  Kawada  and  Mori* 
has  applied  a similar  method  to  droplet  growth  measurements 
using  light  transmission  as  well  as  Mach-Zehnder  Interfer- 
ometry. More  recently  Kalra*  has  investigated  the  condensa- 
tion of  water  vapor  in  an  expansion  wave  using  both  laser 
Fed>ry-Perot  Interferometry  as  well  as  differential  interfer- 
ometry. 

Since  the  flow  history  as  observed  at  a fixed  location 
is  different  depending  on  its  distance  from  the  diaphragm 
location,  the  expansion  wave  offers  the  possibility  of 
studying  the  effects  of  cooling  rates  on  condensation  pro- 
cesses by  making  observations  at  different  locations  along 
the  driver  section  thus  covering  a remge  of  cooling  rates 
in  a single  experiment.  The  average  cooling  rates  in  most 
supersonic  nozzle  condensation  experiments  are  10*  C/s^* 
and  in  Barschdorff 's  shock  tube  experiments  5x10*  C/s.  In 
the  experiments  described  in  this  work  the  typical  average 
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cooling  rates  of  5x10 ^ C/s  were  lower,  approaching  those 
of  the  cloud  ch^unber. 

THE  CENTERED  EXPANSION  FAN 

The  heart  of  the  shock  tube  technique  used  in  the 
condensation  experiments  here  lies  in  the  understanding 
of  the  gasdyncunics  associated  with  the  unsteady  expansion 
flow  in  the  driver  section.  It  is  recalled  that  a simple, 
isentropic,  centered  expansion  wave  is  produced  in  the 
driver  section  provided  that:  (i)  the  complete  removal  of 
the  diaphragm  is  instantaneous  and  the  flow  is  set  into 
motion  abruptly;  (ii)  there  are  no  three-dimensional  and 
viscous  effects;  and  (iii)  the  fluids  employed  behave  like 
perfect  gases.  It  is  well-known  that  in  such  an  isentropic 
expansion  wave  the  flow  properties  of  fluid  velocity,  den- 
sity, pressure,  temperature,  etc.  are  constant  along  the 
characteristics  of  the  governing  differential  equations**’**. 
In  an  X - t plane , these  characteristics  are  straight 
lines  of  constant  x/t,  whose  slope  corresponds  to  the  local 
wave  speed,  c,  given  by 


where  a^  * sound  speed  of  the  undisturbed  medium  in  the 
driver  section, 

Y^  * ratio  of  specific  heats  of  the  undisturbed 
medium  in  the  driver  section, 

p^  s pressure  of  the  undisturbed  medium  in  the 
driver  section, 

p - local  pressure  within  the  expansion  wave  along 
the  characteristic  in  question. 


The  static 
fixed  location. 


pressure  variation  with  time,  p(t),  at  a 
*obs»  given  by. 


P(t_)  _ 

P4 


2Y4 


(2) 


and  the  time  taken  by  the  first  characteristic  (the  head 
of  the  expansion  wave)  to  reach  is  given  by. 


^obs  " *obs^^4 


(3) 


The  isentropic  relationships  from  Poisson's  Equations 
relate  the  remaining  flow  variables  by 


2Y4  2Y4 


Each  fluid  element  traverses  the  expansion  fan  with 
a local  velocity,  u,  given  by 

u s (c  - a.)»  (5) 

^4 

along  a particle  path  given  by 

Y4-I 


*4^ 

p 

‘ - TyvIT 

2 - (Y4-I) 

{^) 

J’ 

(6) 

where  x defines  the  initial  position  of  the  particle  path, 
o 

However,  possible  non-idealities  in  the  real  flow  due 
to  the  bulging  of  the  diaphragm  under  pressure  before 
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rupture,  finite  opening  time,  three-dimensional  initial 
flow,  boundary  layer  effects,  etc.  may  cause  the  real 
expansion  flow  to  deviate  from  that  predicted  by  the 
ideal  theory.  Moreover,  for  condensation  studies,  one 
has  to  trace  each  fluid  element  under  observation  back 
to  its  initial  position  in  the  driver  section  and  inte- 
grate the  eunount  condensed  along  the  particle  path  in 
order  to  compute  the  total  condensate  mass  fraction, 
g at  a given  time  and  at  a given  location, 

Xobs*  order  to  define  the  particle  path  uniquely  the 

starting  point  or  origin  of  the  expansion  fern  has  to  be 
known.  Therefore,  prior  to  the  condensation  study, 
experiments  on  the  flow  itself  using  dry  gases  were  made 
on  the  expansion  flow  actually  produced  in  the  shock 
tube . 

The  experiments  were  performed  in  a 76.2  mm  (3 
inches)  I.D.  steel  shock  tube  which  was  later  also  used 
for  condensation  studies.  Prior  to  an  experiment,  the 
tube  was  outgassed  by  evacuating  it  to  10"*  torr  for  em 
hour  with  a cold-trapped  diffusion  pump.  The  tube  was 
then  filled  with  high-purity  scientific-grade  argon. 

After  recording  the  initial  temperatures  amd  pressures 
in  the  two  section  the  flow  was  started  by  rupturing 
the  stretched  diaphragm  with  a spring-loaded  needle. 

Pressure  time  histories  at  three  different  locations 
situated  at  127  mm  (5  inches),  178  mm  (7  inches),  euid 
495  mm  (19.5  inches)  respectively  upstream  of  the  dia- 
phragm (i.e.  into  the  driver  section)  were  measured  by 
three  Kistler  Model  606L  pressure  transducers  mounted 
flush  with  the  inner  tube  wall  and  electrically  insulated 
from  it.  The  charge  outputs  from  the  transducers  were 
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converted  to  amplified  voltage  signals  by  Kistler 
Dual  Mode  Amplifiers  (Model  504E)  and  were  displayed  on 
the  storage  screen  of  a fast  single  sweep  oscilloscope 
(Tektronix  Model  5103N). 

This  procedure  was  repeated  for  many  different 
values  of  initial  diaphragm  pressure  ratios  ranging  from 
p^/pj^  = 2.0  to  P4/P]^  = 8.0.  Whereas  the  initial  temper- 
atures in  both  sections  were  kept  at  ambient  temperatures, 
the  initial  pressure  in  the  driver  section  was  varied  from 
one  to  two  atmospheres  while  that  in  the  driver  section 
never  went  below  150  torr.  Under  these  conditions,  the 
terminating  characteristics  of  the  relatively  weak  expan- 
sion waves  produced  would  remain  in  the  driver  section. 
Hence  the  outflow  from  the  driver  section  remains  sub- 
sonic. 

Figure  1 shows  a typical  oscillogreun  from  which  points 
of  equal  pressure  at  different  locations  amd  times  were 
mapped  in  an  x - t plane  as  shown  in  Figure  2 . The  iso- 
bars drawn  through  these  experimental  points  are  found  to 
be  straight  lines  merging  at  a single  point  not  far  from 
the  origin  of  the  x - t diagram.  This  result  seen  in 
Figure  2 indicated  that  viscous  effects  were  negligible 
since  a boundary  layer  would  give  rise  to  an  effectively 
variable-area  duct  causing  the  characteristics  to  curve. 
Moreover,  the  interaction  between  a wall  boundary  layer 
and  an  initially  centered  wave  would  lead  to  a non-centered 
wave  ^ • . Figure  3 shows  the  remarkable  agreement  between 
the  wave  speeds  determined  from  the  slopes  of  the  isobars 
and  the  ideal  centered-wave  theory.  This  result  further 
demonstrates  the  fact  that  the  real  expansion  flow  pro- 
duced can  be  viewed  effectively  as  a simple  centered  wave 
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with  a "virtual"  origin  shifted  with  respect  to  the 
actual  diaphragm  location.  The  virtual  origin  can  be 
determined  best  by  noting  the  value  of  which  gives 

a best  fit  between  the  experimental  pressure  variation 
and  the  theoretical  curve  of  Equation  (2).  The  values 
of  t^jjg  for  each  of  the  three  observation  stations  de- 
termined in  this  manner  agreed  well  with  those  found 
from  Figure  2.  Using  experimental  values  of  , the 

pressure  measurements  from  Figure  1 can  finally  be 
represented  by  a single  curve  of  p/p^  versus  as 

shown  in  Figure  4. 


It  is  worthy  of  note  that  the  relatively  low  ini- 
tial diaphragm  pressure  ratios  (2  £ P4/P1  £ employed 
in  these  experiments  produced  only  weak  expansion  waves 
(0.32  ;<  P3/P4  £ 0.70).  Therefore  the  flow  occurred  at 
high  Reynolds  numbers  (Re).  It  is  well-known*^  that 
deviation  from  theory  grows  with  increasing  expansion 
wave  strength,  namely  low  p^  and  low  Re,  which  would  pro- 
duce more  serious  viscous  effects  resulting  in  thicker 
boundary  layers.  Our  primary  error  results  from  the 
diaphragm  that  actually  bulges  under  pressure  before  rup- 
ture, thus  physically  contributing  to  a shift  of  the 
starting  point  of  the  expansion  wave  to  the  right  in 
Figure  2. 


CONDENSATION  EXPERIMENTS 


The  experiments  on  the  condensation  of  H2O  and  D2O 
were  performed  in  the  same  shock  tube  as  shown  schemati- 
cally in  Figure  5.  As  before,  prior  to  an  experiment  the 
whole  tube  was  evacuated  to  10" ‘ torr.  The  condensing 
substance  normally  in  a liquid  state  was  first  injected 
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into  a vaporizing  chamber  heated  by  electrical  heating 
tapes.  The  vapor  was  then  carefully  metered  into  the  shock 
tube  via  a needle  valve  which  was  electrically  heated  to 
avoid  condensation.  The  H2O  used  was  triply  distilled  water 
and  in  the  case  of  D2O,  a 99.8%  Deuterium  heavy  water  was  used. 
(Wilmad  Glass  Co. , Inc.  of  New  Jersey).  The  initial  partial  pres- 

ure  of  the  condensing  vapor  was  measured  with  a calibrated 
Statham  electric  pressure  transducer  while  the  tempera- 
ture was  measured  with  six  calibrated  thermocouples  dis- 
tributed along  the  driver  section.  The  tube  was  next 
filled  with  filtered  and  cold-trapped  ultra-high-purity 
argon  (min.  purity  99.999%).  The  pressure  in  the  driver 
section  also  filled  with  argon  was  preset.  The  gas  mix- 
ture in  the  driver  section  was  then  stirred  for  15  minutes 
by  circulation  through  a sealed  Metal  Bellows  pump. 

After  a steady  state  of  all  properties  was  reached  the 
experiment  was  initiated  by  rupturing  the  diaphragm. 

Two  Kistler  transducers  measured  the  static  pressure 
variation  at  415  mm  and  810  mm  upstream  from  the  diaphragm. 
Rayliegh  light  scattering  at  90° of  a 15  mW  Helium- Cadmium 
laser (A  = 4416  Angstrom)  was  measured  by  two  photomulti- 
plier tubes  (RCA  7265)  having  a fast  response  and  a high 
sensitivity  at  the  wavelength  of  the  laser.  The  pressure 
and  light  scattering  signals  were  recorded  by  an  oscillo- 
scope triggered  by  a third  Kistler  transducer  mounted  50  mm 
upstream  of  the  first  observation  station. 

Onset  of  condensation  was  indicated  by  a pressure 
bump  amd  the  simultaneous  first  appearance  of  the  light 
scattering  signal. 

The  heat  addition  due  to  condensation  in  the  subsonic 
flow  field  generates  waves  which  can  propagate  both  upstream 
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and  downstream,  thus  causing  the  pressv^e  profile  to 
deviate  from  the  isentropic  pressure  variation.  Inside 
the  condensation  zone  the  characteristics  become  curved 
due  to  increasing  sound  speed  caused  by  heat  addition. 

If  condensation  onset,  i.e.  the  critical  supersaturation, 
is  the  major  interest  the  problems  due  to  heat  addition 
can  be  avoided.  The  initial  diaphragm  pressure  ratio 
can  be  adjusted  such  that  the  tail  of  the  expansion  fan 
produced  coincides  with  the  condensaiton  onset  locus  at 
the  observation  station  as  shown  schematically  in  Figure 
6.  Here  the  isentropic  flow  is  preserved  at  all  loca- 
tions in  the  driver  section  up  to  the  point  of  condensa- 
tion onset,  a state  that  coincides  with  the  tail  of  the 
centered  expansion  fan* . Since  the  condensation  onset 
conditions  are  not  known  before  hand  an  iterative  experi- 
mental procedure  is  required  to  "home  in"  on  the  onset 
conditions.  Figure  7 shows  such  an  Onset  experiment 
where  pressure  and  light  scattering  traces  from  an  expan- 
sion with  argon  only  under  the  same  experimental  conditions, 
are  superimposed  on  the  experiment  where  condensation  has 
just  begun.  In  analyzing  the  data  in  such  a situation  it 
is  reasonable  to  assume  that  the  accumulation  of  conden- 
sate and  hence  the  heat  addition  due  to  condensation  are 
negligible  in  the  expansion  fan  prior  to  the  instant  of 
detectable  onset  of  condensation.  Therefore  no  departure 
from  ideal  flow  has  teiken  place. 

RESULTS  AND  DISCUSSION 

Experimental  results  are  presented  in  Table  1 and  con- 
densation onset  conditions  are  presented  in  a plot  of  pres- 
sure versus  temperature  as  shown  in  Figure  8.  Here  the 
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equilibrium  vapor  pressure  curves  of  both  H2O  and  D2O 
as  well  as  some  condensation  onset  data  of  H2O  from 
other  sources  in  supersonic  nozzles**"*’ 

and  a shock  tube*  are  shown. 
Apparently  no  condensation  results  of  D2O  are  to  date 
available  in  the  literature.  In  Figure  8 we  note  a par- 
tial overlap  between  this  work  and  Barschdorf f * s shock 
tube  data.  In  general  the  nozzle  results  lie  to  the 
left  of  the  shock  tube  results.  This  difference  is  to 
be  expected.  It  can  be  explained  on  the  basis  of  a dif- 
ference in  cooling  rates.  The  average  cooling  rates  were 
10*  C/s  in  nozzles,  5xio*  C/s  in  Barschdorff ' s work  and 
5x10*  C/s  in  this  work.  It  is  well-known  that  for  a given 
initial  relative  humidity  greater  supersaturation  can  be 
achieved  at  higher  cooling  rates*®***.  The  same  argument 
also  sheds  light  on  the  observation  (Table  1)  that  the 
onset  data  measured  at  the  first  station  have  a greater 
supercooling  compared  to  those  at  the  second  station.  The 
first  observation  station  is  closer  to  the  diaphragm  lo- 
cation and  therefore  it  has  a cooling  rate  roughly  twice 
that  at  the  second  station  as  can  be  deduced  from  data 
such  as  those  of  Figure  2 . The  described  cooling  rate 
effect  was  found  for  both  H2O  and  D2O  condensation. 

In  order  to  compare  experiments  with  nucleation  theory 
one  needs  to  evaluate  the  condensate  mass  fraction  accumu- 
lated along  a particle  path  finally  reaching  the  experi- 
mentally determined  point  of  condensation  onset . The  for- 
mulation of  this  problem  leads  to  a system  of  integro- 
differential  equations  similar  to  those  first  proposed  by 
Oswatitsch**  for  nozzle  flow.  These  equations  have  been 
treated  elsewhere***®  and  they  will  not  be  repeated  here. 
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The  comparison  between  theory  and  experiment  is  finally 
accomplished  in  terms  of  a single  empirical  adjustment 
factor,  r,  such  that  the  actual  nucleation  in  a given 
experiment,  Jgjjpj  is  given  by^  * 


'^exp  = ^ 


cZ' 


(7) 


where  is  the  nucleation  rate  based  on  the  classical 
theory  of  homogeneous  nucleation  developed  by  Becker  and 
Dttring  (1935)“,  Volmer  (1939)“,  and  Frenkel  (1946)“, 
et  al.  For  each  experiment,  the  value  of  F is  varied  in 
the  corresponding  calculation  until  “ 10”* 

at  the  instant  of  time  when  the  computed  from  Eq.(7) 

matches  well  with  the  measured  onset  conditions.  Since 
there  is  doubt  as  to  whether  the  state  of  the  condensate 
at  onset  conditions  is  of  a solid  or  supercooled  liquid, 
the  experimental  results  were  evaluated  for  both  cases  of 
solid  and  liquid  condensates.^* 


Despite  the  lack  of  reliable  property  data  below  the 
melting  point,  (e.g.  surface  tension  of  solid  H2O  and  D2O), 
and  the  uncertainty  of  extrapolation  when  supercooled  liq- 
uid property  data  are  not  available  (e.g.  D2O),  it  can  be 
seen  from  Table  2 that  the  values  of  logj^QF  obtained  in 
this  work  are  in  reasonable  agreement  with  those  found  in 
previous  studies*’**.  The  relatively  low  values  of  F 
further  confirm  the  remarkable  fact  that  the  classical 
homogeneous  nucleation  theory  in  the  case  of  H20  and  D20 
serves  as  a good  prediction  of  the  onset  of  condensation. 
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FIGURE  CAPTIONS 


Figure  1 


Figure  2 
Figure  3 
Figure  4 

Figure  5 
Figure  6 

Figure  7 


Oscilloscope  traces  showing  experimentally 
measured  pressure  variations  at  three  different 
locations  situated  at  127  mm,  178  mm  and  495  mm 
respectively  upstream  from  the  diaphragm  loca- 
tion. For  the  first  (bottom)  and  second  (top) 
stations,  99  Torr/div.  vertical;  for  the  third 
station  (middle),  109  Torr/div.  vertical.  The 
common  time  base  is  0.5  ms/div.  horizontal. 

The  gas  employed  was  argon,  p^  = 777.1  Torr; 
P4/P1  = 3.70;  T^  = Tj^  = 297.4  K. 

Pressure  measurements  of  Fig.  1 mapped  into  an 
x-t  plane. 

Wave  speeds  determined  from  the  x-t  diagram  in 
Fig.  2.  Centered  Expansion  Wave  Theory. 

p/p4  versus  plot.  O first  station; 

□ second  station; ^ third  station;  centered 

expansion  wave  theory. 

Schematic  diagram  of  the  shock  tube  for  conden- 
sation experiments.  (Dimensions  in  mm.) 

Condensation  zone  within  the  centered  expansion 
fan.  Pressure  variations  at  *nd  along  a 

particle  path  reaching  the  point  of  condensation 
onset. 

Oscilloscope  traces  showing  isentropic  flow  in 
the  expansion  fan  is  preserved  in  a condensation 
experiment  up  to  the  point  of  condensation  onset. 
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Figure  8 

Table  1. 
Table  2. 


Plot  of  pressure  versus  temperature  showing 
experimentally  determined  condensation  onset 
conditions.  # H2O,  this  work;  ■ D2O,  this 
work;  ^ H2O,  Nozzles  ‘ ; Q H2O  shock  tubes* 


Initial  and  onset  conditions  of  a selected 
number  of  experiments . 

Summary  of  H2O  and  D2O  condensation  onset  data. 
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Table  2.  Summary  of  H2O  and  D2O  condensation  onset  data 


Substance 

AT=T. 

K mp 

log 

10  ^ 

References 

(°C) 

solid 

liquid 

Steam 

0 to  + 40 

-2  to  1 

Nozzles  (Experi- 
ments**”** and  r 
calculations*  * ) 

H2O 

-62  to  -19 

5 

to  7 

3 to  4 

Nozzles  (Experi- 
ments^”"** and  r 
calculations*  * ) 

H2O 

-48  to  -12 

— 

1 to  6 

Shock  Tube* 

H20 

-41  to  -16 

1 

to  3 

3 to  6 

This  work 

D2O 

-51  to  -21 

1 

to  2 

1 to  3 

This  work 

Figure  1 
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Condensation  of  Water  in  Argon 


Superposed  with  Expansion  in  Pure  Argon 
Pressure 

T 

155Torr/cliv 
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